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We propose an entanglement classification for symmetric quantum states based on their diagonal matrix- 
product-state (MPS) representation. The proposed classification, which preserves the stochastic local operation 
assisted with classical communication (SLOCC) criterion, relates entanglement families to the interaction length 
of Hamiltonians. In this manner, we establish a connection between entanglement classification and condensed 
matter models from a quantum information perspective. Moreover, we introduce a scalable nesting property 
for the proposed entanglement classification, in which the families for N parties carry over to the N + 1 case. 
Finally, using techniques from algebraic geometry, we prove that the minimal nontrivial interaction length n for 
any symmetric state is bounded by n < \_N/2\ + 1. 


Entanglement is widely considered the cornerstone of 
quantum information and an essential resource for relevant 
quantum effects, such as quantum teleportation [1-4], quan¬ 
tum cryptography [5, 6], or the speed-up of quantum com¬ 
puting [7], as in Shor’s algorithm [8]. Moreover, entangle¬ 
ment is recognized as useful for understanding properties of 
condensed matter models, such as quantum phases [9] and 
topological orders [10], among others. Entanglement based 
properties are usually challenging to study both experimen¬ 
tally and theoretically, due to the exponential growth of the 
associated quantum degrees of freedom. Experimentally, due 
to the exponentially large amount of degrees of freedom typi¬ 
cally involved, the detection and the quantification of the en¬ 
tanglement are difficult to achieve. Theoretically, quantities 
describing the entanglement are generally complicated func¬ 
tion of the quantum state, and they are normally arduous to 
analyse. To overcome these obstacles, advanced quantum in¬ 
formation techniques have been successfully applied to an¬ 
swer condensed matter questions [11-13], shedding a distinct 
light on the problem. This novel approach may bring about 
exciting results in many-body physics, and result in a new 
revolution in physics, where quantum information and matter 
phenomena can be formally unified [14]. 

An important question in quantum information is the clas¬ 
sification of entanglement by means of some mathematical or 
physical equivalence. Classifying entanglement should help 
in recognizing similarity between different entangled states, 
and it may be useful to boost the practicabilities of quantum 
information protocols. A first result is that quantum states 
connected by SLOCC operations, which perform probabilis¬ 
tically the same quantum tasks, can be collected into entan¬ 
glement classes, called SLOCC classes, but also known as 
SLOCC criterion [15]. Nevertheless, there is an infinite num¬ 
ber of SLOCC classes for four or more parties that may be 
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gathered, in turn, into a finite number of entanglement fami¬ 
lies [16-21]. Unfortunately, the community has not been able 
to relate all classes and families to specific properties or quan¬ 
tum information tasks, although a few of them have certainly 
raised experimental interest [22-25]. It is noteworthy to men¬ 
tion that, up to now, no general characterization nor classifi¬ 
cation of entanglement exist for many-body systems. 

In this Article, we present an entanglement classification 
for quantum states induced by their MPS structure, which pre¬ 
serves the SLOCC criterion and is exemplified for the sym¬ 
metric subspace. The proposed classification is based on the 
local properties of the multipartite quantum state. In this 
sense, it relates entanglement families to the interaction length 
of Hamiltonians, establishing a connection between entangle¬ 
ment classification and condensed matter. Our proposal is 
twofold beneficial: on the one hand, it does not result in an 
infinite number of entanglement classes, considerably simpli¬ 
fying their study; on the other hand, it provides a direct phys¬ 
ical insight to the nonlocality of entanglement classes, given 
by the interaction length of their parent Hamiltonians. Addi¬ 
tionally, we introduce a scalable nesting property in which the 
families for N parties carry over to the W + 1 case. 

We focus on the classification of the SLOCC classes corre¬ 
sponding to symmetric states, which are invariant under any 
permutation of the parties [27], i.e. Fltf) = \tjj) where F is 
an exchange operator. This is an interesting subspace, since 
its dimension grows linearly with the number of parties but, at 
the same time, it contains many physically relevant states. 

A pure state 14/) is called entangled if it is not separa¬ 
ble [26], i.e. if it cannot be written as a tensor product 
14/) |4/i) ® 14 / 2 ) ® • • • 0 |4/Ar). This definition is indi¬ 

rect and this may be the reason that its quantification for more 
than two parties is still unsettled. The idea of separability 
emerges when one tries to identify which states can be gen¬ 
erated from other states, defined locally in each subsystem, 
by using local operations and classical communication among 
patties. Therefore, these local operations provide a natural cri¬ 
terion to collect quantum states in entanglement classes with 
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the same type of entanglement [15]. More specifically, two 
states belong to the same class if they can be transformed into 
each other with non-zero probability via SLOCC operations. 

The splitting of the Hilbert space into the SLOCC classes 
is fully understood for bipartite and tripartite qubit systems, 
even in the nonsymmetric case [15]. There is only one entan¬ 
glement class for two qubits; the one containing the Einstein- 
Podolski-Rosen (EPR) state |T') = ^ (|00) -f |11)), and 
two symmetric classes for three qubits: one represented 
by the Greenberger-Horne-Zeilinger (GHZ) state |T') = 
^ (|000) + I 111)), and the second one referenced to the W 

state |T') = ^ (|100) -I- |010) -f |001)). However, for four 
or more qubits, the number of SLOCC classes is infinite, and 
their parametrization grows exponentially with the number 
of parties, while lacking robustness against experimental er¬ 
rors [15]. In this sense, SLOCC classification makes the asso¬ 
ciation of classes to specific physical properties difficult. This 
explains, so far, the lack of experimental interest in producing 
states beyond the well-known GHZ or W, among few others. 

Due to the natural interest of the SLOCC criterion, it is 
customary to collect these infinite SLOCC classes into a fi¬ 
nite number of larger families based on specific mathematical 
properties [16-18, 20] or operational approaches [21, 28] (see 
Fig. la). However, up to now, all classifications have failed 
to associate groups of classes or families to specific quantum 
information tasks. Since there are infinite SLOCC classes for 
four or more parties, they can be partitioned into families in 
an infinite number of ways, and we would not expect all to 
be relevant. To overcome this conundrum, we propose the 
following criteria that an SLOCC classification into families 
should fulfill: 1) Every SLOCC class must belong to only one 
family (classes should not cross any border separating fami¬ 
lies), 2) separable states must be in one family and this fam¬ 
ily should contain only separable states, 3) SLOCC classes 
belonging to the same family must show common physical 
and/or mathematical properties, and 4) the classification into 
families must be efficient in the sense that (a) the number of 
families must grow “slowly” with the number of qubits, and 
(b) the efforts for classifying N qubits should be useful for 
classifying -f 1 (nesting). 

In the last decades, matrix product states (MPS) and 
tensor networks have emerged as a powerful tool to 
tackle complex problems in many-body systems [29- 
31, 33]. In this sense, any quantum state |T') = 

Sii iiv=i • ■ ■’^n) can be rewritten in a local man¬ 
ner as IT-) = • • ■ 4!^’ 1*1 • ^n), where 

are matrices containing all local information related to site k. 
This language is convenient to describe ground states of lo¬ 
cal Hamiltonians [30-32], sequential processes [34-36], and 
systems fulfilling an area law [37]. In the following, we use 
the MPS formalism to define an entanglement classification of 
quantum states into families of SLOCC classes, based on the 
local dimension of the matrices describing the states, called 
bond dimension (see Fig. lb). Furthermore, we prove that this 
MPS classification is directly related to the interaction length 
of the corresponding parent Hamiltonian [31]. We apply these 
novel concepts to the case of the symmetric subspace, al¬ 


though they could be extended to more restricted or more gen¬ 
eral sets of states. We start by highlighting that any symmetric 
state admits MPS representations with site-independent diag¬ 
onal matrices. If we further request a minimal bond dimen¬ 
sion, we must find the optimal way to represent any symmetric 
state as The dimension of the ma¬ 

trix parallels the number of the Schmidt coefficients, which 
has been proposed to quantify entanglement [38]. Hereafter, 
when we refer to the bond dimension of a symmetric state, we 
mean the minimal bond dimension associated with its diago¬ 
nal representation. 

Crucially, SLOCC transformations preserve the minimal 
bond dimension (see Supplementary Information). Indeed, if 
IT'^^) and are two quantum states with bond dimen¬ 

sions Da and Db respectively, and C a class, then 

= (1) 

This implies that all states of the same SLOCC class may be 
represented with the same minimal matrix dimension, which 
is a SLOCC invariant. In this way, we can define a family of 
SLOCC classes by means of the following equivalence rela¬ 
tion; 

Definition (Diagonal MPS entanglement classification). Let 
Sa and Sb be SLOCC classes, and Da, Db the minimal bond 
dimension of their respective states in its diagonal MPS rep¬ 
resentation. We say that Sa ~ Sb Da = Db, and we call 
entanglement families the resulting equivalence classes. 

Notice that the class of separable states can be optimally 
represented with matrices with bond dimension D = 1, and, 
indeed, it coincides with the family D = 1. Therefore, the 
diagonal MPS (DMPS) classification, proposed here for the 
symmetric subspace, fulfills the first two aforementioned con¬ 
ditions. Moreover, in this DMPS classification any symmetric 
state of N qubits can be written with at most bond dimension 
N, so the number of families grows linearly with the number 
of parties involved. In this sense, the DMPS classification also 
satisfies a recently proposed tractability criteria [20]. 

The explicit translational invariance of the MPS 
formulation leads the DMPS classification to ful¬ 
fill the aforementioned criterion (4b). Let {Ai}i 
define an iV-partite symmetric state |\I/(*^)) = 

Ei'str |ii...ijv). Then, the state 

which lives in a different Hilbert space, namely that of 
iV -f 1 parties, does show exactly the same local properties 
as |T'(^)). As the DMPS classification respects these local 
properties, and not just for a given state but for the whole 
SLOCC class, an intriguing nesting property of the families 
for different number of parties emerges; 

Theorem (Nesting). Let us consider an N-particle symmetric 
state of qubits j^’Ar) = \xk)^^ with optimal bond 

dimension D^fjj^), such that D^ippj) < \_N/2\ + 1. Then, 

the state |'0jv+i) = has optimal bond 

dimension D{'il)]si+i) = D(jPn). 
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FIG. 1: a, The SLOCC criterion divides the Hilbert space (the square) in such a way that every quantum state is in a well defined class 
(the lines). For four or more parties, the number of these SLOCC classes is infinite. However, they may be gathered into families (the 
colored areas) under certain rules, ideally with more physical associations than mathematical ones. Here, the condition is given by the 
minimal bond dimension of the matrix-product-state representation of quantum states, relating the MPS classes to the interaction length of 
parent Hamiltonians, b. The proposed MPS classification enjoys a scalable nesting property in which the classes of an A'^-partite family can be 
mapped onto the classes of the corresponding {N -f 1) case, generating a malryoushka structure. A detailed example is given for the symmetric 
subspace of arbitrary number of parties. 


See Supplementary Information for the proof. This theo¬ 
rem shows that, from the perspective of the local properties, 
the only purely TV-partite states are the ones whose optimal 
bond dimension is larger than the maximal bond dimension of 
any state with TV — 1 parties. This generates a matryoushka 
structure depicted in Fig. lb, where, unlike other entangle¬ 
ment classifications, the classification for TV parties is con¬ 
nected with the classification for TV + 1, for all TV. We be¬ 
lieve that a further exploitation of this scalable nesting prop¬ 
erty would be interesting in the many-body case, where the 
exact number of particles is usually not relevant. 

Lastly, in the MPS formalism, the role of parent Hamiltoni¬ 
ans for a given state comes to the fore. Namely, for any given 
state, one can construct a local Hamiltonian which includes 
it in its ground eigenspace. The MPS representation informs 
us about the features and interaction length of this construc¬ 
tion. For instance, it controls whether it is the only ground 
state or there is a spontaneous symmetry breaking [29], or the 
inheritance of local and global symmetries [30, 31]. This con¬ 
structive method works as follows: Let |vl/(^)^ be a symmet¬ 
ric pure state of TV parties with an MPS representation with 
bond dimension D. The reduced density matrices for n < N 
are defined by = tvN-n i.e. by trac¬ 

ing out N — n parties. By construction, rank(p^”)) < D, 
so when n > log 2 D, has a kernel. Let > 0 be 
the projector onto this kernel and H = the to¬ 

tal Hamiltonian. Thus, iT is a positive operator and |vl/(tv)^ 
is in its ground manifold since = 0. Additionally, 

the interaction length of H is n. In order to apply this con¬ 
struction to the DMPS classification, notice that all reduced 
density matrices for more than one party have nontrivial ker¬ 


nel when acting on the full n-qubit Hilbert space, since their 
support is restricted to symmetric states. We must then con¬ 
sider the relevant kernel in the corresponding symmetric sub¬ 
space of dimension n -F 1. Clearly, if the bond dimension is 
D, the rank of is at most D in the symmetric space as 
well. Thus, we have ensured the existence of a relevant par¬ 
ent Hamiltonian with interaction length n, if n > D. How¬ 
ever, techniques from algebraic geometry allow us to prove 
that the minimal nontrivial interaction length for any symmet¬ 
ric state is bounded as n < [TV/2J + 1 (see Supplementary 
Information). We shall see from the examples below, indeed, 
that interaction lengths shorter than D do arise. Notice that 
in our construction we propose families of parent Hamiltoni¬ 
ans, for all of which all the states in an entanglement family 
are ground states. In particular, the ground manifold is always 
degenerate in our construction. To put this statement in con¬ 
text, bear in mind that the concept of parent Hamiltonian of a 
state is a general one: a local Hamiltonian which includes the 
state in its ground manifold . A particularly well established 
technique for constructing parent Hamiltonians for states or 
families of states is given by the MPS construction. In that 
case, the uniqueness of the ground state of this MPS Hamil¬ 
tonian will depend on certain properties of the corresponding 
matrix. In particular on the property (or otherwise) of injectiv¬ 
ity, such that a lower bound for the interaction length is given 
by a polynomial in D and d [39]. A particular example of this 
for the GHZ state may be found in Ref. [40] 


i) GHZ states .— The GHZ states may be represented, in¬ 
dependently of TV, by the matrices Aq = diag{l,0} and 
Ai = diag{0,1}. The corresponding parent Hamiltonian 
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family reads 


N N 

Hqhz = J ^ ^ CTi ■ <^i+l ~ Jz ^ '^i+l t 

i=l i=l 

with the conditions > 0 and > 2 J. 

ii) W states .— The Wn states generalize to N par¬ 
ties the structure of the three-partite W state, and 
they correspond to Dicke states with just one excita¬ 
tion. They can be represented with bond dimension 
D{Wn) = N hy the following sequence of diagonal 
matrices Aq = A diagjO, 1, ..., ..., 

with (fi = 2 tt/N{N — 1), and Ai = pdiag{(l — 
1, ..., ..., where a = 

2ttIN. Then, we propose the parent Hamiltonian family 

N 

Hw=JY1 - T'T'i ■ 0■^+l] , 

i=l 


where both J and 7 need to be positive. Notice that, even 
though the W state requires a bond dimension N to be repre¬ 
sented with diagonal matrices, we can find a Hamiltonian with 
interaction length 2 which has this state as a ground state. 

Hi) X]sf states .— Let us consider the family for TV >4 
given by 




(iV-1) + z^-^Vn\Wj^) 


up to normalization. The DMPS representation of this 
state is given by the {N — 1)-dimensional matrices Aq = 
z diag { 1 , ..., e 2 ’^AJV- 2 )/(Ar-i) | 

Ai = diag{1,1,1,...!}. In this case, the relevant interac¬ 
tion length is seen to be 3. In fact, the parent Hamiltonian 
family reads 

Hx = ^ J2Z1 (3 + Serf crf+iCT|+2 - «+i - erf crf+2 
-crf+iCrf+2 - Ej=o <+i) + i EZi (1 - ■ ^^+l) , 

for any J > 0 and > 0. 

In conclusion, the proposed MPS classification fulhils all 
the aforementioned criteria for a versatile grouping of SLOCC 
classes, while maintaining a linear growth in the number of 
families with respect to the number of parties. At the same 
time, the unveiled nesting property allows us to use, in a scal¬ 
able manner, the effort invested in the MPS classification for 
N parties in the subsequent TV -f 1 case. Additionally, we have 
provided a physical meaning for the proposed MPS classes by 
connecting them to paradigmatic properties in condensed mat¬ 
ter. This missing link in the theory of entanglement classifica¬ 
tion has been exemplihed for the case of the symmetric sub¬ 
space, where the DMPS representation was used. It is note¬ 
worthy to mention that the provided DMPS classification can 
be formally extended from the symmetric subspace to more 
general sets of quantum states. We believe that MPS-based 
entanglement classifications will be able to connect mathe¬ 
matical aspects already known in quantum information theory 
with relevant physical features in many-body systems of ex¬ 
perimental interest, such as quantum state preparation or the 
emergence of fractional magnetizations. 
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Appendix A: The symmetric subspace and SLOCC classes 

We consider the symmetric subspace of TV-particle systems, 
i.e. Sym (7/®^). Eor the sake of clarity, we say that a state 
is symmetric if and only if F\il}) = |' 0 ), where F is a flip 
or exchange operator, or the representation of any element of 
the permutation group. If the one-particle Hilbert space % has 
dimension d (i.e. qudits), we have 

In particular, the dimension of the symmetric subspace in an 
TV-qubit system is TV4-1. Our goal is to organise in some phys¬ 
ically relevant way the orbits of the group of invertible local 
operator (ILO) transformations, i.e. the orbits under transfor¬ 
mations of the form 


IV'tv) , (AI) 

where A is an invertible local operator acting on H. The 
reason for considering orbits under the group of ILOs is that 
stochastic local operators and classical communication oper¬ 
ators (SLOCC) are constructed as ILOs [15]. Notice that, in 
the symmetric case, we only need to consider ILOs of the form 
above, i.e. such that the invertible local operator acting on one 
party is the same for all parties [41]. 
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Appendix B: Diagonal MPS representation 


Proof. Let us write the one-site states in the standard basis; 


The main result of the paper relies on the diagonal matrix 
product state (DMPS) representation of a symmetric state. In 
this section, we give two independent proofs that any symmet¬ 
ric state can be written in this form. 


1. First proof: Majorana representation 

Our starting point is the Majorana representation: for N- 
qubit systems any symmetric state can be written as 

IV'tv) CX ^ Isii) ® lOa) ® ® |ei„) , (Bl) 

l<il^i2¥=-"¥=iN<N 

where 

|e,) = a,|0)+/3,|l) (B2) 

are single qubit states. In the case of qubits, this decomposi¬ 
tion is unique up to rearrangement of indices. Then, Eq. (Bl) 
can be rewritten as 


\xk) = I3k\l) + ak\0) (B 6 ) 

with \akf + |/3feP = 1. At most one of the /J^’s can be zero, 
since otherwise we would have the state | 1 ) repeated, and the 
set would not be pairwise linearly independent. 

On the other side of the equivalence, if more than one of the 
/?fe’s are zero, the set of one-particle states are not 

pairwise linearly independent. Let us assume first that none 
of the /3fe’s are zero. 

The A^-site state can be written as 

\xk)®^ = (/3fe|l)+afe|0))®'^ 

Linear independence of the set means 

that the only set of -|- 1 coefficients Xk such that 
= 0 is Afc = 0 yk. However, any linear 
superpositions of can be written in terms of 

the Dicke states as 


/N-k 

ctGPjv /c=0 ^ ^ \ /=! / 

This state is of the form 

D 

IV'Af) CX \xk)'^^ , (B4) 

k^l 

for some D, that, as we will see, can be easily cast in a diago¬ 
nal MPS representation. 


iV+l N 

y] = (B8) 

k—1 n—0 

with the following relation between the { 7 ^} and the {Ai} co¬ 
efficients; 

7» = \/(f) (E j . (B9) 

This means 


2. Second proof: Vandermonde determinant 

We give now an alternative proof of the existence of diago¬ 
nal presentations (B4) for qubits. The dimension of the sym¬ 
metric subspace is W -|- 1 , thus, if we can show that A^ -|- 1 lin¬ 
early independent vectors of the form \xk)^^ exist, we have 
shown that any state can be written as their linear combina¬ 
tion. We will find useful in what follows to relate states to 
their expansion in terms of Dicke states: 

N 

= (B5) 

k=0 


/7o\ 

71 


X \ 

A2 

\inJ 

= SVB 

VAat+i/ 


(BIO) 


where S = diag ( 1, 


diag (/3f, /3f,..., and 





where these qubit Dicke states are defined by the number of 
“excitations” k, i.e., in the standard basis {|0), |1)}, the num¬ 
ber of times 11) appears, and they are normalised. 

First, we prove the following statement; 

Lemma 1. The set of N-particle states is lin¬ 

early independent and generates the whole symmetric space 
if and only if the one-particle states are pairwise 

linearly independent one particle states. 


V 1 1 ••• 1 / 

The condition for linear independence of the set 
is thus equivalent to the condition for V 
to be invertible. Notice that ^ is a Vandermonde matrix, and 
its determinant is 


1^1 


n 

i<fe<i<Ar-l-i 


ai. 

Pi 


ak 

Pk\ ’ 


(B12) 
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which is different from zero if and only if the set 
is pairwise linearly independent. We have thus proven the 
lemma in the case that none of the jSk coefficients are zero. 

Let us now consider the possibility that one of the /? is zero, 
e.g. /?! = 0. Then, we have w.l.o.g. ai — 1, and 


/7o\ 

7i 

= SVB 

a‘ ] 

\jn/ 


VAat+i,/ 


(B13) 


where B = diag (l, 13^,..., and 



(B14) 


Linear independence of the set at hand is equivalent to the 
invertibility of V, which reads 


/ \ ^-1 

/ \ ^-1 


' ttiV + l \ 

[P2) 

V ^N+l J 

1 

1 

= n2<fc<i<Ar+l 

Oil o^k 


Pi Pk 



(B15) 


Since we have required that all j3 be different from zero, save 
for Pi, we recover again the requirement that the one site vec¬ 
tors be pairwise linearly independent, thus com¬ 
pleting the proof. □ 


Corollary 1. IfM -f 1 > D, then the set of M-particle states 
{\xk)^^'\^^i is linearly independent if and only if the one- 
particle states {\xk)}k^i are pairwise linearly independent. 

It follows that any symmetric state can be written in the 
form of Eq. (B4) with D < N -\- 1, where the vectors 
{|®fe>}£i^ are pairwise linearly independent. This gives us 
also a bound for the optimal D, i.e. for the minimal D in the 
decomposition of a state. In what follows, we use D{'f) to 
denote the optimal D of the state \ more than one state is 
being discussed. Otherwise, we normally denote the optimal 
value by D without further qualification. Thus, as a corollary 
of Lemma 1 we conclude that for all N qubit symmetric states 
D <N + 1. 


3. Diagonal Matrix Product States (DMPS) 

In the two previous subsections, we have shown that a rep¬ 
resentation of the form given by Eq. (B4) is always possible 
for any symmetric state of N qubits. In this subsection, we 
show that this is equivalent to a diagonal MBS representation. 


Any quantum state |<1>) S admits an MBS representa¬ 
tion given by 

2 

|$)= ^ Tr[AW.-.4^1]|/ri.../r^) (B16) 

where the set of complex matrices 1C = € 

= 1)2, 1 < fc < N} are called Kraus operators 
or MPS matrices and D is the bond dimension of the MBS 
[42]. 

If the state is translational invariant, then there always exists 

\k] 

a site-independent MBS representation such that = Ai 
Vfc, normally by increasing the bond dimension of the MBS 
[29][Theorem 3]. Moreover, two sets of Kraus operators {Ai} 
and {Bi} are equivalent [29] [Theorem 2], i.e. they represent 
the same state, if there exists an invertible matrix X such that 
Bi = XAiX-^,yi. 

It is clear that MBS representations implemented by diago¬ 
nal matrices lead to symmetric states. That the opposite also 
holds is suggested in [29][Apendix A.l]. We state it as a The¬ 
orem, whose proof is immediate given the previous construc¬ 
tions. 

Theorem 1. Any N-qubit permutation invariant state 
= J2k=i ak\xk)®^ has a diagonal MPS representa¬ 
tion, with bond dimension D < N -\- 1, given by the Kraus 
operators 


D 

Af, = '^{ck)^ {p\xk) \k){k\, 

k=l 


(B17) 


where D is the optimal bond dimension if and only if the de¬ 
composition = 'Ylik=i^k\xk)®^ uses also the mini¬ 

mal number of vectors. 

Proof. Let us replace the Kraus operators given by (B17) in 
(B16): 

2 

= E \pl---PN){pl---PN\Xki---Xk„)^ 

ki,...,kN 

A. A. 

^^ki * ‘ ‘ ^k]^^kik2^k2k^ • ■ • ^k^ki 

= ^Ck\Xk)®^- 
k^l 

The relation between the optimal bond dimension and the 
minimal number of vectors in the representation is trivial by 
construction. This concludes the proof □ 

After this existence check, it now behoves us to ask for the 
optimal bond dimension for a given state is, and after its prop¬ 
erties. 
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Appendix C: Entanglement and optimal bond dimension 

In this section, we justify the dehnition of entanglement 
family by proposing and discussing the basic properties that 
we require. Moreover, we compare this dehnition with those 
in the literature, and we prove that our classihcation is stable 
with respect to changes in the number of parties. 

1. Definition of family and comparison with previous works 

W. Diir, G. Vidal and J. I. Cirac established in 2000 that the 
relevant notion for entanglement classes is equivalence under 
SLOCC [15]. Moreover, they proved that SLOCC transfor¬ 
mations are implemented by ILOs. Indeed, |^) and \(j>) are 
connected via an SLOCC if and only if a set of local invert¬ 
ible operators exists, such that 

N \ 

IV')- (Cl) 

i=i / 

If we are restricted to the symmetric subspace, the existence of 
symmetric ILOs connecting two states jV') and \(j)) is enough 
to ascertain their SLOCC equivalence; i.e. they are SLOCC 
equivalent if there exists a local invertible operator A such that 
[41]. 

One should notice that Local Unitary (LU) equivalence is 
included in SLOCC equivalence. Thus, any LU entanglement 
classihcation is fully included in one and just one SLOCC en¬ 
tanglement class. In this sense, and according to the dehnition 
of “family” that follows, we can say that SLOCC entangle¬ 
ment classes are actually/ami/iei with respect to LU classes. 

Definition 1 (Entanglement family). An entanglement family 
in a set of states is a collection of entanglement classes for the 
same set with the following properties: 

• An entanglement class fully belongs to one and only one 
family (SLOCC structure preservation) 

• The class of separable states must be in a family of its 
own. 

• The number of families grows in a controllable manner 
with the number of particles. 

• There is a mathematical criterion for the arrangement 
of classes in families, with a physical interpretation. 

Optionally, one might add an additional desideratum: 

• Independence of the criterion with respect to the num¬ 
ber of parties N. 

As stated, the SLOCC classification can then be understood 
as an arrangement in entanglement families of the LU entan¬ 
glement classes. Even so, in what follows we will refer to 
SLOCC classes, because they are more relevant in what re¬ 
gards to quantum information tasks. In the following, we 
propose our sectioning into families and compare it to other 
previous classifications. 


a. DMPS classification 

We propose that the optimal D diagonal bond dimension be 
the criterion for classification. It is well known, from the form 
(B4), that it is an SLOCC invariant, and it has been proposed 
as a measure of entanglement [38]. Here, we examine its value 
as a classification tool. 

As D is an SLOCC invariant, an entanglement SLOCC 
class belongs to one family, and just one. All separable states 
have D — 1, and if a state presents D — 1 then it is separa¬ 
ble, by definition. Thus the separable class lies in a family of 
its own. As we see, the first two criteria we have proposed 
for families are met. In the next sections, we will discuss 
the physical interpretation of our classification, and its depen¬ 
dence on N. 


b. Comparison to the degeneracy classification 

The organisation in families proposed in T. Bastin et a/ [21] 
relies on the Majorana representation. We assert that this clas¬ 
sification in families meets our requirements. Eor all N the 
family Vm includes the separable class and only the separa¬ 
ble class. The degeneracy configuration which is the classifi¬ 
cation criterion is an SLOCC invariant. The third property is 
elucidated in L. Lamata et al [43]. Instead, the fourth tentative 
desideratum is not so clearly fulfilled. 

Notice that this classification is different from ours. In fact, 
already for N = 4 their family includes classes that 

would go to our D = 2 and D = 3 families. 

c. Comparison to SLOCC normal form classification 

In E. Verstraete et al [44], they give a procedure for con¬ 
structing normal forms under SLOCC operations. It is a clas¬ 
sification of SLOCC classes into families, since all the ele¬ 
ments of a SLOCC class will go into the same normal form. 
It does not, however, meet all our desiderata: the states in the 
Wn class have normal form equal to zero, which is also the 
normal form for separable states. Thus, the class of separable 
states is not in a family of its own. 

d. Comparison to entanglement polytopes classication 

In Ref. [20], Walter et al provide an organization of SLOCC 
classes in families according to polytopes which classify the 
one particle reduced density matrix. Eor the specific case of 
symmetric states of N qubits, which they compute, they ob¬ 
tain [A^/2] -f 1 families; i.e. growing linearly with the number 
of qubits, thus fullfilling our requirement of controlled growth 
of the number of families. Separable states however belong 
to all entanglement polytopes (although this might be taken as 
quibbling, since they are the only states for which the relevant 
parameter is 1). An operational characterization is given by 
the maximal linear entropy of entanglement achievable in a 
family. 
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2. Nesting and stability 

In this subsection, we prove a limited Theorem on nesting, 
in order to clarify how our classification is stable with respect 
to changes of N. Indeed, we define the stability property as 
follows: the family is stable if, when moving up in number of 
parties, one maintains the minimal bond dimension and thus 
the local properties. To clarify: since we consider the family 
for N qubits defined in terms of the optimal bond dimension 
of its MPS representations, we move up in number of parties 
by starting with an A^-qubit state in the family, fixing its ma¬ 
trices Aq and Ai, and now using the selfsame matrices for 

+ 1 qubits. This can be done for all states in the family. We 
prove now that stability is a generic property, in that it holds 
for D < [N/2\ + 1, given N. We shall explain later (Sub¬ 
section D 3) in what sense this is indeed generic. Since we are 
dealing with minimal bond dimensions for different states (in 
different Hilbert spaces, even), we explicitly denote the opti¬ 
mal bond dimension for state |'0) by D{ij}) in this context. 

Theorem 2 (Nesting). Let us consider an N-particle sym¬ 
metric state of qubits I^Pm) = with optimal 

bond dimension such that D('ip]\j) < [A^/2J -f 1. 

Then, the state \ipN-\-i) = has optimal 

bond dimension = Z?('0jv). 

Let us introduce the Schmidt binary rank at level M for the 
studied systems, which we will denote, for an A^-particle sym¬ 
metric state |'0Ar), as sbrj^f (V'Af)- It is defined as the Schmidt 
number for the bipartition of the system in M and N — M 
qubits, i.e. as the minimal number s such that \iIjn) can be 
written in the form 

S 

iv-w) = \£,k)M\Xk)N-M , (C2) 

fc=l 

with {|Cfc)M}fc=i orthonormal vectors in the space of M 
qubits, and similarly {|xfc) N-M}k=i in th® space of N — M 
qubits. The following Lemmas are useful to prove the Nesting 
Theorem. 

Lemma 2. sbr^ {iPn) < Dlppjq) for all 1 < M < N. 

Proof We can write [iPn) = Yhk^"'^ \xk)^^■ By defin¬ 
ing |6 )m = and |xfc)Ar-M = \xk)’^^^~^\ and 

comparing with Eq. (C2), it directly follows that, for all 1 < 
M < N and all symmetric W-particle states, sbrM (V’w) ^ 
DiipN) ■ □ 

Lemma 3. If \iPn) = Yl,k=i \ik)M\Xk) n-m, where 
{\ik)M}k=l and {|Xfe)A''-M}£=i o.re two sets of linearly in¬ 
dependent vectors, then shTMif^N) = D. 

Proof There exist two invertible operators Xa.b, such that 
the sets {|ffc)M}£=i and {\Xk)N-M}k=i’ where = 

XA\^k)M and \xk)M = XB\xk)N-M, are orthogonal. 

This implies that sbrM('0Af) = D for the state IV'tv) = 
Xa 0 XB\f}N)- As the Schmidt binary rank is invariant un¬ 
der local invertible operations, we have that sbrM(V'Af) = 
sbrM(^/’Ar) = D. □ 


Proof of Theorem!. We have that D(fjM+i) < D{ijj]s[) by 
construction. On the other hand, if DfipM) < then 

an integer M exists such that Diyipjsi) < M-\-l and Diytpjsi) < 
N -\- I — M. As {\xk)}^^’^'’ are necessarily pairwise lin¬ 
early independent. Corollary 1 implies that, in the bipartition 
\^N+i) = (|xfe)®“) 0 (|xfe)®(^+t-^)), both sets 

and {linearly 
independent.Then, Lemma 3 ensures that sbrM (V'at+i) = 
Lemma 2 now implies that D{f)A!) < O('0jv-i-i)> 
which concludes the proof. □ 


3. Examples of DMPS representation 

In this subsection, we provide examples of optimal matrix 
product representations. This comes with a general method to 
find the optimal bond dimension for a given state. 


a. Statement of the problem 

Let us consider an W-particle symmetric state expressed in 
the Dicke basis 


N 

= (C3) 

Q—0 

as our initial data. In order to identify a Matrix Product 
representation, we shall use the equivalence of Theorem 1. 
Thus, we look for D one-qubit states of the form \xk) = 
a:fc|0) -I- yk\f), such that we can write \tj)) = E^=i \xk)^^ ■ 
This demands 


T.^k-^y'^k = {\ d^, (C4) 

k=l ^ d 

for all a from 0 to N. The task at hand is to examine, given 
the state, which is the minimum D for which the set of equa¬ 
tions (C4) has a solution, and, once having identified such 
minimum D, to obtain such solution. 

It is important to notice that, even if we were to start by 
positing two D x D diagonal matrices Aq and Ai, which do 
determine a symmetric N qubit state, we would not have de¬ 
termined whether that bond dimension is indeed optimal for 
the state. By construction, the optimal value would be smaller 
than or equal to the bond dimension proposed, and it is still 
pertinent to examine the set of solutions of Eq. (C4). 

In the following examples, we give optimal D construc¬ 
tions, that are clearly not unique for a given state. 


b. Wn state 


The Wn state is defined as 


\Wn) = ^(|100... 000) -f |010 ... 000) + • • • -f 
-f |000...010)-f I000...001)). 
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with N > 3. Equivalently, this is the Dicke state with one in Eq. (C4) are thus da = 5ai, with 5ap the Kronecker delta, 
excitation, namely The coefficients da to be inserted Let us write a solution in terms of matrices; 


ylo = Adiag { 1 , e 2 WtV(W-l)^ . . . ^ ^ 2 .^k/N(N-l) ^ _ ^ 2 .^(N- 2 )/N(N-^) 

= Mdiag {l, e 2 -dtv-l)/tv^ , , , ^ ^ 2 MN-k)/N^ ^ ^ (1 _ ^)l/tv| _ 


(C5) 

(C6) 


This is adequate since, for 0 < Z < iV, we have 


Tr [A^-^A\] 


N-2 


= E 


^27Tik{N-l)/N{N-l)^277il{N-k)/N 


fc=0 


N-2 

_ ^2mi ^ ^2Trik{l — 1)/{N — 1)] 

k=0 

1 — exp [27rj(l — 1)] 

1 — exp [27ri(l — 1)/{N — 1)] 

= <5a(ZV-l), (C7) 

and Tr [^i^] = 0. Therefore, we can choose any A, /r 
respecting the normalisation relation A^(/r/A)(7V — 1) = 
1/y/N. Additionally, we can also fix the gauge, by means 
of Al|Ai cx 1. One can now write this solution in terms of 
\xk) states using Eq. (B17) with Ck = 1: 

\xk) = /r|l)+ , l<fc<iV-l 

\xn) = 

The existence of this form guarantees that D{Wn) < N. We 
need to prove that indeed D{Wn) = N. In order to do that, 
consider the set of equations (C4), particularised for this case, 
for a generic D. Assume first that none of the coefficients 
Xa is zero. Then, up to normalisation, the equations can be 
written in the equivalent form 


1 1 

2:1 ^2 




1 ' 

ZD 




/<\ 


N 


/ \<j 


0 

\oJ 


(C8) 


where Za = yajxa- Let us define the (N + 1) x 79 matrix 


1 1 


where Za are all distinct, and the + 1 component vector 


h = 


/o\ 

1 

0 

Voy 


(CIO) 


The requirement that the rank of Zd be equal to the rank of 
the augmented matrix [Zd b] is a necessary condition for the 
existence of solutions to the set of linear equations (C8) for 
Xa, with Xa 7 ^ 0 Vo!. As rank (Zd) — 79 (because all the 
Za must be distinct), it is a necessary condition for existence 
of solutions of the system (C4) that the rank of the augmented 
matrix [Zd b] also be 79. 

The case in which one of the coefficients Xa = 0 can be 
written in a similar manner. Without loss of generality, as¬ 
sume that the null coefficient is the 79-th one. Then, the sys¬ 
tem of equations can be written as 


/ 1 1 


1 \ 


-21 


^2 

Z 2 


ZD-1 

^D-1 




/ X^ \ 




-ij 


r^N 

J^2 


\^D-lJ 


= b- 


0 

Vz/^y 


= bD- 


(Cll) 

Then, a similar condition is retrieved, namely that the rank of 
Zd -1 and the rank of the augmented matrix are equal, which 
is a necessary condition for the existence of solutions for the 
system of equations. 

Let us start by examining the case of 79 = TV — 1, initially 
for Xa =/= 0 Va. The necessary condition for the existence of 
solutions is that all the TV x TV minors of 

/ 1 •• • 1 0 \ 

Zi • • • Zn-1 1 

0 


^N-l 


w 


(C12) 


-1 oy 


are zero, using the hypothesis that for all a ^ /3 it is the case 
that Za ^ zp. Thus, compute hrst 



Zl 

Z 2 • 

■ Zd 


-Zl • 

■ ZN-1 

1 



1 • 

■ 1 

Zd = 

4 

z 2 . 

Z-2 

■ 4 

, (C9) 

z 2 . 
Zl 

■ 72 

^AT-l 

0 

= (- 1)^+1 

(n^) 

2:1 • 

• Zn-1 



Z^ • 

Z 2 

'■ 4 


Z^ ■ 
Zl 

• Z^ 

ZN-1 

0 



N-2 

Zl 

N-2 

ZN-1 
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Since the last factor is a Vandermonde determinant, in order 
for this minor to be zero, it is necessary that one of the Za is 
zero. Without loss of generality, let that be = 0. We 

then compute 

N-X N-\ 

Z-x ^N-2 

1 • • • 1 

ZN-2 


In order for this minor to be zero as well, we would need an¬ 
other Za to be zero. That is however excluded, since we re¬ 
quire from the beginning that all the to be different. It 
follows that the two minors computed so far cannot be simul¬ 
taneously zero, while keeping pairwise linear independence of 
the {Since all iV x At minors must be zero for the 
rank of the augmented matrix to be At — 1, and two of those 
cannot be simultaneously zero, under the assumption that all 
Xa 7 ^ 0, it follows that the rank of the augmented matrix is 
larger than the rank of Zn-1 and there is no acceptable solu¬ 
tion of this form. 

Let us now relax the assumption that all Xa 7 ^ 0, and al¬ 
low xn-x = 0. Then, a similar analysis with the augmented 
matrix [Zn -2 bN-i], where b is defined in Eq. (Cll), shows 
again the inexistence of solutions. 

The arguments above hold for any D < N—1, as the image 
of Zjji is included in the image of Zjy, if D' < D. We have 
thus finally proven that the optimal bond dimension for Wn 
is indeed N. 

By means of this example, we have provided an algo¬ 
rithm for the determination of the optimal bond dimension 
for a symmetric state of N qubits: 1) construct the b vec¬ 
tor from the Dicke basis expansion of the state, and define 
bfy = b — {0,0, ■ ■ ■ , 2) start with £> = 1; 3) Is it possi¬ 

ble for rank[Z^ b] to be D or for rank[Z^_j^ 6^] to be Z? — 1? 
If yes, then D = D. If no, 4) let D be set to D + 1 and return 
to 3. 


n 

,Q=1 


2:1 


M-3 


1 ••• 1 10 

Zl • • • ZN-2 0 1 

Zx ■■■ z %_2 0 0 

• • • z^~^ 0 0 

1 ^N-2 


c. GHZm state 


It is obvious by construction that the GHZat states 
IGHZat) = (lOO---O)-f 111-••l))/y2 have D = 2 and do 
not really require the machinery presented above. Nonethe¬ 
less, it can be instructive to apply it to this case. We have to 
distinguish between two cases: a) A^ = 2 and b) A^ > 2. Let 
us first tackle A^ = 2. In this case we have to study whether it 
is possible for 



1 

-Z2 


(C13) 


to be of rank 2. This would require that the determinant is 
zero, i.e. {z 2 — zx) (1 + Z 1 Z 2 ) = 0 , which, under the con¬ 
dition Zx 7 ^ Z 2 entails Z 2 = —fjzx- Inserting this into the 
system of equations, we have 

x\zrx\ = 1, (C14) 

zxxf -^ = 0, (C15) 

^1 

z!xf + ^ = 1. (C16) 

The second equation provides us with = zfxf. Substitut¬ 

ing in the others we are led to a;f = 1/2 and zf = 1, which in 
turn give us a;| = 1/2. This is a valid solution, as xx ,2 7 ^ 0. 
We recover thus the following D = 2 decomposition of the 
GHZ 2 state: 

IGHZ 2 ) = ^ (|0) + |1))®^ + ^ (|0) - | 1 ))®^ . 

(C17) 

Passing now to the general case, W > 3, we first attempt 
D = 2 with the hypothesis xx,X 2 7 ^ 0. This entails examining 
the rank of 


/I 1 IX 
Zx Z 2 0 
zf 0 

W 4 1/ 


(CIS) 


We should study all 3 x 3 minors, in principle. It is however 
enough to look at just two: first 


1 

Zx 


1 1 
Z2 0 
0 


ZXZ 2 {Z 2 - Zx) , 


(C19) 


which, without loss of generality, leads us to require Z 2 = 0 , 
since we require it to be zero, and that zx 7 ^ Z 2 - Under this 
Z 2 = 0 condition we next examine 


1 

-21 


1 1 
0 0 
0 1 


-Zi . 


(C20) 


For this to be zero as well we would have to require zx = 0, 
which we cannot allow by our requirement that zx 7 ^ 22 - Thus 
there is no H = 2 presentation for IGHZat) with N > 3 such 
that both xx and X 2 are different from zero. Let us now assume 
X 2 = 0. We have to study the system 





Zx 

N 

0 


Xi = 






0 

U/ 


The solution is given by zx = 0, Xx = 1, 
corresponds to the trivial construction 


(C21) 


1, which 


|ghz^) 


_(|0)®A^ + |l)®tV 


(C22) 
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d. Xn state 


Xa 7 ^ 0 Va. We need that all the {N — 1) x {N — 1) minors 
of 


We now give an example of a family of symmetric N- 
particle states {N > 4) that is explicitly in a different SLOCC 
class with respect the previous ones. Let us denote it by 
Xm{z). They are given hy D = N — \ one-qubit vectors 
as follows: 


/ 1 ••• 1 0 \ 

N—l 

Zi • • • Zn-2 z' 

7^ ... 7^ 0 

^iV-2 


(C25) 


\xu) = |1) + , (C23) 


Vi 


z 


N 

N-2 


1 


with k = 1,..., N — 1 and z 0. The state is thus, by a 
non-normalised representative, 

1 w-i 

^ k^l 

= \l)^^ + z^-^s/N\Wn) ■ (C24) 

The case = 4 (with z = 2“^/®) was introduced in A. 
Osterloh and J. Sieweit [45] as a maximal entangled genuine 
four-partite state, with the notation |$ 2 ). These authors also 
introduced the general case N, with z = [46]. 

In the following, we prove that the optimal bond dimension 
for Xff{z) is N — 1, by applying the same method as in the 
Wn and GHZ^v cases. The case D = N —1 has a solution by 
construction. Let us examine the case of D = N — 2, first for 


are zero, in order to have solution. Let us compute first 


2^1 • 

• Zn-2 

zN-1 




z! ■ 

■ Z%_2 

0 

1 

1 — 1 

II 

/N- 

■ n 

zf- . 

■ z^-1 
ZN-2 

0 








1 ••• 

1 





Zi 

ZN-2 




X 

^N-3 

Zl 

^N-3 

ZN-2 


zi X 


. (C26) 


In order to be zero, it is necessary that one of the Zq, be zero. 
Let us assume, without loss of generality, that z^_2=0. We 
now compute 


1 

■ 1 

1 

0 






1 • 

■Zl 

■ 1 

• Zm-3 

^1 

• 

• Zn-3 
■ Z^ 

Zn-3 

0 

0 

ztv-i 

0 

= z^-i 

. 

• z2 

Zn-3 

= z^-i 

(ff^) 






N-2 

N-2 


\a=l / 








Zl 

■ ^tV-3 


N-A 

N-A 

N-2 

Zl 

N-2 

Zn-3 

0 

1 





Zl 

■ zjv_3 


(C27) 


In order to be zero, we would need another Zq, to be zero, but 
this case is excluded by hypothesis. If we assume now that 
X]si-i = 0, then a similar argument brings us to the conclusion 
that the optimal bond dimension is W — 1. 

Appendix D: Parent Hamiltonians 
1. Introduction 

Given a state ip we say that a Hamiltonian H is a parent 
Hamiltonian of ip if this state is a possible ground state of H, 
not necessarily unique [31]. For finite-dimensional systems, 
a trivial parent Hamiltonian can be built by simply consider¬ 
ing the projector onto the state ip, denoted by P^. In fact, for 
any positive a, the Hamiltonian H = a {1 — Ptf,) has ip as its 
only ground state. The concept of parent Hamiltonian comes 
in as useful when it can be written in terms of short range or 
very well-controlled long-range interactions. In the MPS rep¬ 
resentation, the considered state is written as locally as it is 


possible, and the minimal bond dimension for an MPS pro¬ 
vides us with some control over the range of the Hamiltonian 
in question. One systematic way of searching for a given par¬ 
ent Hamiltonian is to consider the M -particle reduced density 
matrices, that, for some number M on, will not be full rank. 
It is enough to identify the projectors onto the kernels, and 
compute the intersection of their complements. The projector 
on the complement to that intersection is a parent Hamiltonian 
for the state at hand, with shorter interaction range. 


2. Parent Hamiltonians for symmetric states 

The case of symmetric states is particularly simple, since 
all M-particle reduced density matrices are identical. 

Let us give the name Pm to the M-particle projector onto 
that kernel of Notice that Pm is an M particle opera¬ 
tor. Acting by translation on Pm 0 (i.e. producing 

1 0 Pm 0 ^ ^ 

cesively), and adding the translated terms, we obtain a parent 
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Hamiltonian for the state at hand, with interaction length M. 

Let denote the reduced density matrix for M qubits 
in this symmetric case, that is = Tin-m {\'>Pn){'4’n\)- 
Assume that IV'w) has optimal bond dimension D, i.e., it can 
be written as IV'at) = J2k=i for some non-normalised 

vectors which are pairwise linearly independent. 

We have that 

p(M) ^ ^ ({Xj\xk))^~^ {\xk){Xj\)’^^ . (Dl) 

It follows that rank < D, since the image of p^^^ is 

included in the span of { \xk)’^^ } ■ This gives us a first bound 
on the interaction range. In fact, acts on the space of M 
qubits, which is a linear space of dimension 2^. If the rank is 
smaller than the dimension, the kernel is not trivial. Clearly, 
if M > log 2 D it follows that the kernel of p^^^ is not trivial, 
and thus we know that the interaction length needed will be at 
most [log 2 D\ +1. 

Notice that the previous bound has been obtained allowing 
to project onto the whole linear space of M qubits, being the 
vectors symmetric or not. If we allow the projection onto the 
orthogonal part of the symmetric space, we can obtain a par¬ 
ent Hamiltonian whose interaction range is two, and whose 
ground state space includes the whole symmetric space. For 
instance, it is enough to consider the projection onto the state 
101) — 110). However, our goal is to build parent Hamiltonians 
able to discriminate amongst symmetric states. Therefore, we 
derive a bound on the interaction length for non-trivial parent 
Hamiltonians, whose ground state space contains the orthog¬ 
onal part of the symmetric space. 


It follows that the rank of is precisely that of the lin¬ 

ear transformation T'm- Furthermore, we have also shown 
that rank = rank[p(^)], since, from equation 

(D3), we can conclude that, as matrices, 

'I'Af-M = 'I'm ■ (D5) 


The interaction length is the longest when the ranks of all 
the reduced density matrices achieve their maximal possible 
value. For values of M from 1 to a maximum to be deter¬ 
mined, the reduced density matrix is full-rank in the symmet¬ 
ric space when rank [p^^^] = M -\- 1. On the other hand, 
since [piN-M)^ = rank [p(^)] in all cases, we have that 
has rank at most 2, at most 3, and so on. Thus, 

the last possible value of M for which it is possible for 
to be full rank is \_N/2\. Namely, 


rank 


^(YN/2\+l) 


= rank 


^(W2l-i) 


< \N/2\ < lN/2\ -f 1. 


(D6) 


Therefore, with certainty, ker |^p(L^/2J+i)j is not empty, and 
we conclude that, for all states, there will be parent Hamilto¬ 
nians with interaction lengths smaller or equal to [A^/2J + 1. 

Notice that the states whose optimal bond dimension is 
larger than [A^/2J will show interaction lengths smaller than 
the optimal bond dimension. Although we shall not expand on 
the concept here, it is true that the set of states with optimal 
bond dimension smaller or equal to [A^/2J -f 1 are dense in 
the space of states, and we may call it the set of generic states. 


4. Examples of parent Hamiltonians 


3. Generic states and sharper bound for interaction lengths 

In the following, we prove the bound n < [A^/2J -f 1 for 
the interaction length of parent Hamiltonians built considering 
only the symmetric space. We prove it by using a technique 
translated from algebraic geometry [47]. Given an A^-qubit 
symmetric state \iIjn) we define for each M from 1 to W — 1 
the linear transformation m that maps M-qubit symmetric 
bra vectors, to {N — M)-qubit symmetric ket vectors 

by 

M {{4>m\) = {4>m\'<Pn) ■ (D2) 

Given bases in the space of symmetric states of M qubits, 
and of N - M qubits, {\e%_M)}k=f^^^ the 
matrix elements for the linear transformation 4/^ are com¬ 
puted to be 

(^M\u={{e'N-M\®{eM\)\^N). (D3) 

The reduced density matrix for N — M qubits is now written 
in the basis as 

= ('f'M)ife fe) . 

A / km 


In all successive examples, we will have to supplement the 
Hamiltonians with a two-particle term that implements the re¬ 
striction to the symmetric space; i.e. the two-particle projector 
onto the spin 0 sector, the singlet. This is given by 

(Po).,.+i = (D7) 

a. GHZm state 

The M-particle reduced density matrix for the |GHZjv) 
state is 

p(^)(GHZAr) = i (|0)(0|®^ -f |1)(1|®^) , (D8) 

which is always of rank two. The one-particle reduced density 
matrix is maximally mixed, while the two-particle reduced 
density matrix has the states |10) and |01) in its kernel. The 
corresponding projector is 

(Ps=l,m=o)i,,+i = ^ (1 + cr, • cr,+i - 2crfCTf+i) (D9) 
= ^ (1 + «+i + • 


(N-M) 

Plm 


(D4) 
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Therefore, we identify the GHZ]v state as a ground state of 
the Hamiltonian 


N N 

-ffCHZjv = > (DIO) 

i=l i=l 

with the conditions > 0 and > 2 J. 

By construction, it is easy to see that this Hamiltonian 
presents degeneracy, since |0)®^ and |1)®^ are both eigen¬ 
states with the same energy. 


For this expression to be valid for M = 1,2, define IIT 2 ) = 
(|01) -f |10)) /y/2 and |VFi) = |1). Therefore, a possible par¬ 
ent Hamiltonian is 


N 


iFi 


w 


Ct V 

~ R 2-^ 


-2crf 


2/3 


O’-,- 

a 


(D12) 


with positive a and j3. Also in this case, the ground state is 
degenerate. 


b. Wn state 


We have 

pW [Wm) = ^|VFm)(1^m| + . (Dll) 


c. Xn state 


The reduced density matrix for the |X 7 v( 2 :)) state in the 
M < N — 1 particles case reads 


[X^{z)] = ^ ^ |^|L-2jv [M\Wm){Wm\ + {N-M) |0)(0|®“] } . (D13) 


Therefore, for all 3 < M < — 1 the rank of the reduced 

matrix is 3. The case M = TV — 1 is different, and it reads 




+ -f . (D14) 


The case M = 2 can be included in the general formula for In order to compute an alternative form for this local Hamil- 
arbitrary M with the convention that tonian term, we look at the following operator: 


|1T2) = ^(|01) + |10)). 


(D15) 


Again, this means that (Xn) is of rank three, which is ac¬ 
tually the size of the symmetric space for two qubits. Thus, the 
only element of the kernel of the two-particle density matrix 
is the projector out of the symmetric space, i.e. the projec¬ 
tor onto the singlet. In order to have a local Hamiltonian that 
discriminates this state, we need to go to the three-particle 
density matrix. The Dicke state with two excitations is the 
only subspace of the symmetric space that lies in the kernel of 
p*^^^ (Xn), and provides us with the local Hamiltonian term 
/i(3) = \d!^^){D^^\, of interaction length 3. 


n 


1 = 0 ,l^k \i=l 


aj — n + 21 


In the symmetric space, it is zero except when acting on 

= 0, on which it has the value HlLo i^k ~ 2fc). 
Thus, we have the projector 


p(") _ 

^k ~ 


n 

l=0,l=^k 


T,i=i < - ^ 
2{l - k) 


(D16) 


which, when restricted to the symmetric space, provides us 
with the relevant vector. On explicit computation for the case 
at hand, we find the relevant operator 


( 3 ) _ 1 


pw; _ 

-^2 ~ 


O ^z I O^z^z^z \ ^ ^2 

3 (72(73 -f 3fT]^ (72(73 / . <7^ 
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which provides us with the parent Hamiltonian presented in 
the main text. 
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